Shell-model half-lives for r-process waiting point nuclei including first-forbidden 

contributions 
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We have performed large-scale shell-model calculations of the half-lives and neutron-branching 
probabilities of the r-process waiting point nuclei at the magic neutron numbers TV = 50, 82, 
and 126. The calculations include contributions from allowed Gamow- Teller and first-forbidden 
transitions. We find good agreement with the measured half-lives for the TV — 50 nuclei with charge 
numbers Z = 28-32 and for the TV = 82 nuclei 129 Ag and 130 Cd. The contribution of forbidden 
transitions reduce the half-lives of the TV = 126 waiting point nuclei significantly, while they have 
only a small effect on the half-lives of the TV = 50 and 82 r-process nuclei. 

PACS numbers: 21.60.Cs, 23.40.-s, 26.30.Hj 



I. INTRODUCTION 

Although the actual site of the astrophysical r-process 
is still not known with certainty, it is commonly accepted 
that occurs in an explosive environment of relatively high 
temperatures (T ss 10 9 K) and very high neutron densi- 
ties (> 10 20 cm" 3 ) 

Under such conditions, neutron 
captures are much faster than competing beta decays and 
the r-process path in the nuclear chart proceeds through a 
chain of extremely neutron rich nuclei with relatively low 
and approximately constant neutron separation energies 
(S n < 3 MeV). Due to the relatively stronger binding of 
nuclei with magic neutron numbers, the neutron separa- 
tion energies show discontinuities at the magic numbers 
TV = 50, 82, and 126. As a consequence the r-process 
matter flow slows down when it reaches these neutron- 
magic nuclei and has to wait for several beta decays 
(which are also longer than for other nuclei on the r- 
process path) to occur until further neutron captures are 
possible carrying the mass flow to heavier nuclei. Thus 
matter is accumulated at these r-process waiting points 
associated with the neutron numbers TV = 50, 82, and 
126 leading to the well-known peaks in the observed re- 
process abundance distribution. 

The beta half-lives of the waiting points have at least 
two important effects on the r-process dynamics and 
abundance distributions. At first, they mainly determine 
the time it takes the mass flow within the r-process to 
transmute seed nuclei to heavy nuclei in the third peak 
around A ~ 200. Second, in the astrophysical environ- 
ment the nuclear r-process timescale (given by the sum of 
beta half-lives of nuclei in the r-process path) competes 
with some dynamical timescale of the environment, e.g. 
the expansion timescale of the ejected matter. If the r- 
process path and half-lives were known, the reproduction 
of the abundance distribution can be used to constrain 



the conditions of the astrophysical environment. If the 
r-process has sufficient time for beta-flow equilibrium to 
establish, the relative elemental abundances are propor- 
tional to the beta half- life [6j. 

Despite their importance only a few half-lives of wait- 
ing points with magic neutron numbers TV = 50 and 82 
are known experimentally while no experimental 

data exist yet for the TV = 126 waiting points. The situ- 
ation is expected to improve in the near future with the 
advent of new experimental facilities. For example, the 
beta-decay half-lives of 38 new neutron-rich isotopes from 
Kr to Tc close to the r-process path have been measured 
at the new RIBF facility at RIKEN [ll|. Furthermore, 
researchers at GSI have measured half-lives of nuclei close 
to TV = 126 using a novel analysis method [l2j. Despite 
this progress, the half-lives needed for r-process simula- 
tions have mainly to rely on theoretical estimates. As the 
Q- values involved are rather low, such calculations have 
traditionally been based on allowed, i.e. Gamow- Teller, 
transitions. Most of these studies used the Quasiparti- 
cle Random Phase Approximation (QRPA) either on top 
of semi-empirical global models [13l - [l5j ] or the Hartree- 
Fock-Bogoliubov method [16[ . Although the calculations 
give a fair account of the few experimental half-lives, it is 
well known that these models underestimate the correla- 
tions among nucleons which pull down the Gamow- Teller 
(GT) strength to low energies. This shortcoming is over- 
come within the interacting shell model which indeed de- 
scribes the measured half-lives of r-process waiting point 
nuclei very well [BUI. 

It is expected that the appearance of intruder single 
particle states with different parity may have influence 
on the low-energy spectra of the r-process waiting point 
nuclei. Thus it is conceivable that first-forbidden tran- 
sitions might contribute to the half-lives of these nuclei. 
A first attempt to estimate such forbidden contributions 
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has been taken within the gross theory [14| . This model, 
however, has been found as rather inaccurate when ap- 
plied to Gamow- Teller transitions. More recently, Borzov 
extended the QRPA studies based on the Fayans energy 
functional to a consistent treatment of allowed and first- 
forbidden contributions to r- process half- lives [13] ■ While 
these calculations find that forbidden contributions give 
only a small correction to the half-lives of the N — 50 
and N = 82 waiting point nuclei, they result in a signif- 
icant reduction of the N = 126 half-lives. This impor- 
tant finding has been our motivation to extend our shell 
model calculations of waiting point half-lives to include 
also first-forbidden transitions. We expect that correla- 
tions among nucleons will not only affect the half-lives, 
but a reliable description of the detailed allowed and for- 
bidden strength function is needed to estimate the prob- 
abilities for beta-delayed neutron emission rates which 
are known to be important to describe the decay of the 
r-process nuclei towards stability after freeze-out. 

We note that GT and higher multipole transitions are 
relevant to describe neutrino-nucleus reactions which are 
important in many astrophysical sites pp . l22j . Tra- 
ditionally, these reactions have been studied within 
the Random Phase Approximation [2^, [24j], includ- 
ing nuclei relevant to r-process nucleosynthesis [25U27j | . 
In an interesting recent development, neutrino- induced 
reactions on li ght nuclei with relevance to neutrino- 
nucleosynthesis [H, Hl| have been calculated on the basis 
of the shell- model, including GT and first-forbidden tran- 
sitions [30l. l3lj. 



II. SHELL MODEL AND ft DECAY THEORY 

In our half-life calculations, we consider allowed and 
first-forbidden contributions. These are obtained using 
the diagonalization shell mode code NATHAN developed 
by Etienne Caurier [H, HH to calculate the initial and 
final nuclear states and the corresponding nuclear tran- 
sition matrix elements. Model spaces and residual inter- 
actions will be discussed for the three different sets of 
waiting point nuclei with N — 50, 82 and 126 individu- 
ally. The partial half-life, t, for a transition between an 
initial (normally the ground state) and a final nuclear 
state is related to the phase space factor by: 



ft = K = 6146 s. 
The phase factor has the form: 



(1) 



C(W)F(Z, W)(W 2 - 1) 1/2 W{W Q - W) 2 dW . 

(2) 

C(W) is the so called shape factor that depends on the 
electron energy, W in units of the electron mass. Wo is 
the maximum electron energy, also in electron mass units, 
that is given by the difference in nuclear masses between 
the initial and final nuclear states, Wq = Q/(m e c 2 ) = 
(Mi - M f )/m e . F(Z,W) is the Fermi function that 



corrects the phase space integral for the Coulomb dis- 
tortion of the electron wave function near the nucleus. 
The partial decay rate is related to the partial half-life: 
A = ln2/i. The total decay rate is given by summing 
over the partial decay rates to all possible final states. 

For allowed transitions, the shape factor does not de- 
pend on the electron energy and for (3~ decay has the 
form: 

C{W) = B(GT). (3) 
The GT reduced transition probability is given by: 
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where the matrix element is reduced with respect to the 
spin operator a only (Racah convention (34{) and the sum 
runs over all nucleons. For the isospin lowering operator, 
we use the convention i_n = p. Finally, (<m/<?v) — 
— 1.2701(25) is the ratio of weak axial and vector coupling 
constants. 

For first-forbidden (FF) transitions, the shape factor 

is: 



C(W) = k + kaW + kb/W + kcW 2 



(5) 



where the coefficients k, ka, kb 1 kc, depend on the FF nu- 
clear matrix elements, the maximum electron energy, Wq, 
and the quantity £ = aZ/(2R) with R the radius of a uni- 
formly charged sphere approximating the nuclear charge 
distribution [35[ . Following the treatment of Behrens and 
Buhring [36[ they are given by: 
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with 



Y = ?y- + O, Ci = Y + |(u - x)W . 



(7) 
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The numbers in parenthesis on the closing bracket de- 
note the rank of the operators inside the bracket. The 
parameter 71 is given by y/l — (aZ) 2 . For the Coulomb 
functions fii and A2 we use the approximations /ii ~ 1 
and A 2 ~ 1 H3. 

After a non-relativistic reduction, the matrix elements 
can be related to the form- factor coefficients, a ' v Fku, 
defined in refs. (36|. l37il . In the Condon and Shortley 
phase convention [38] the matrix elements are: 
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where 



C'lm — \l 21+1 ^ m ' 



(9) 



with Yim the spherical harmonics. The weak axial cou- 
pling constant is $a — —1.2701(25) and M is the nu- 



cleon mass. The quantity 7(1, 1, 1, l,r) appearing in 
the primed matrix elements takes in account the nuclear 
charge distribution that can be approximated by a uni- 
form spherical distribution (36[: 
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(10) 

Based on the conserved vector current theory and the 
assumption that the isospin is a good quantum number, 
the matrix element can be related to the x matrix 
element 13911: 



gy = E 1 x, 



(11) 



where energy E 1 is defined as the energy difference be- 
tween the isobaric analog of initial state and the final 
state: 



E-y = E ias(l) -E f = Q + AE C - (m n c 2 - m v c 2 ), (12) 

where m n and m p are the neutron and proton masses 
and AEq is the Coulomb displacement energy between 
isobaric analog states that can be approximated by pol ]: 



AE C = 1.4136(1) Z/A 1/3 - 0.91338(11) MeV, (13) 

with Z = (Zi + Z f )/2. 

To compare the first-forbidden and Gamow- Teller 
transitions, we define the averaged shape factor: 



C(W) = f/f , (14) 
where / takes the form of Eq. [2] and fo is 

fo = j F(Z,W)(W 2 -l) 1/2 W(W -W) 2 dW. (15) 



III. MODEL SPACES AND QUENCHING 

We have performed beta-decay half-live calculations 
for r-process waiting points based on large-scale shell 
model calculations. In particular we chose the following 
model spaces and respective interactions. 

For the N = 50 nuclei we have adopted a model space 
spanned by the 0/7/2,5/2 and 1^3/2,1/2 orbits for protons 
and by the 0/ 5 / 2 , 1^3/2,1/21 an d 0g g / 2 orbits for neutrons. 
The single-particle energies and the residual interaction 
are the ones adopted in [4l[ to study the shell evolution 
between 68 Ni and 78 Ni. 

Our shell model calculations for the N = 82 waiting 
point nuclei follows the shell model studies presented in 
[l9( . From the two model spaces defined in |19j we adopt 
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the one built on a 88 Sr core. That is we explicitly consider 
the lpi/2 proton orbit which is expected to be impor- 
tant for the description of the negative parity states and 
hence the first-forbidden transitions. Our model space 
is then spanned by the 0g 7 / 2 , 1^3/2,5/2. 2s 1/2 , 0h n/2 or- 
bits outside the N = 50 core for neutrons, and the 
lPi/2j0g 9 /2,7/2) 1^3/2,5/2, 2s 1/2 orbits for protons. This 
model space avoids spurious center-of-mass excitations 
by omitting the 0hn/2 orbit for protons and the 0gg/ 2 
orbit for neutrons. We adopt the residual interaction 
given in [l9| based on 88 Sr core which gives a good ac- 
count of the spectroscopy of nuclei in the neighborhood 
of 132 Sn. In particular, our calculation reproduces the 
excitation energy of the first 1 + state in 130 In as well as 
the low-energy spectrum of 128 Cd and of the r-process 
waiting point nucleus 130 Cd [42j . 

The model space for the N — 126 waiting points has 
been spanned by the 0g7/ 2 , 1^5/2.3/27 0^n/2 and 2si/ 2 or- 
bits for protons and the 0/i 9 / 2 , 1/7/2,5/2, 0ii 3 / 2 , 2p 3 / 2 ,i/ 2 
orbits for neutrons. As interaction we use the effective 
Kuo-Herling interaction KHH e of ref. [43[ which has been 
constructed based on holes in a 208 Pb core. It is the 
same model space and effective interaction as has been 
used in a previous calculation of the half-lives, which, 
however, has only considered pure GT transitions [2lj |. 
These are mainly connected to neutron 0hg/ 2 to proton 
O^-n/2 transitions. However, it is expected that first- 
forbidden transitions can compete, mainly via neutron 
0*13/2 to proton 0hn/ 2 transitions. Full diagonalization 
in this model space exceeds current computer capabili- 
ties. Hence, we performed truncated calculations follow- 
ing a generalize seniority scheme that allows for config- 
urations with maximum seniority 8, i.e. we consider a 
maximum of 4 non-J = pairs, for even even nuclei. 
For odd-even nuclei, the number of broken pairs had to 
be limited to three, while for 199 Ta no limitation has 
been enforced. We expect our model spaces to be large 
enough to give a reasonable account for the low-lying 
Gamow- Teller and first-forbidden transitions. Neverthe- 
less the model spaces are too restricted to recover the full 
Gamow- Teller and first-forbidden strengths built on the 
ground or isomeric states. The missed strength, however, 
resides mainly outside the Qp window and does hence not 
affect our half-life calculations. 

Although the shell model usually gives a good account 
of the relative strength distributions, it overestimates 
the total strength. For Gamow- Teller transitions this 
shortcoming can be corrected for by replacing the bare 
Gamow- Teller operator by an effective operator GT e g- = 
q at. The quenching factor q has been found to be ap- 
proximately constant over the nuclear chart [H]-[46j]. In 
practice, using q sa 0.7 has been shown to give a good re- 
production of the absolute Gamow- Teller distributions. 
There is evidence that also the absolute first-forbidden 
transition strength is overestimated within shell model 
approaches. Ejiri and collaborators [13, EH related this 
fact to core polarization effects and suggested the intro- 
duction of a constant hindrance factor. Based on per- 



turbation theory, Warburton [43|, Ell showed that the 
quenching of first-forbidden transitions depends slightly 
on the initial and final single particle orbits. In partic- 
ular, Warburton found that transitions mediated by the 
rank operator (Eq. 8h) (called relativistic matrix ele- 
ment) appear to be enhanced compared to the other first- 
forbidden transitions, due to meson exchange effects [Hoj . 
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FIG. 1. Comparison of calculated first-forbidden average 
shape factors, obtained for the best fit values of the quenching 
factors (Eq. (16)), with experimental data [5ll45^ ]. 

To treat the quenching of the Gamow- Teller and first- 
forbidden transitions in our shell model calculations we 
have assumed that the quenching factors are the same for 
all nuclei. Following the findings of Warburton [4!|[5(| we 
have furthermore assumed that the quenching factors for 
the operators of rank 0,1, and 2 contributing to the first- 
forbidden transitions as defined in Eq. (8) can be differ- 
ent. To determine these individual quenching factors we 
have performed shell model calculations for experimen- 
tally known beta-decays of nuclei in the vicinity of the 
magic neutron numbers N = 82 and N = 126. Here, we 
have adopted the similar set of first-forbidden transitions 
of nuclei in the lead region as chosen in the study of War- 
burton (50[, supplemented by the decays of the ground 
state of 205 Au (N = 126) and the (1/2) - isomeric states 
in 131 In (N = 82) and 129 In (N = 80) which are both 
known to decay by first-forbidden transitions. By per- 
forming a least-squares fit to the experimental data we 
obtained the following quenching factors for the various 
matrix elements defined in Eq. (8): 

q(£'v) = 1.266, q(w) = q(w') = 0.66, 

q\x) = q{x') = 0.51, q(u) = q[u') = 0.38, (16) 

q\z) = 0.42. 

The calculated half-lives and the corresponding average 
shape factors are summarized in Table[l] Fig. [T]comparcs 
the experimental and calculated shape factors. 

Upon a closer inspection there is quite a good agree- 
ment between our shell model half-lives for the isomeric 
states in the In isotopes with the experimental values. 
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TABLE I. Comparison of calculated log f a t and (c(W) 
for first-forbidden transitions with experimental data [5ll45 



1/2 




Transition 



log fot 



C(W) 



1/2 



Initial 


Final 


the 


exp 


the 


exp 


131 In(r) 


131 Sn(| + ) 


5.32 


«5.1 


65.8 


85.4 




131 Sn(i + ) 


5.74 


6.5 


41.1 


17.1 


129 In(H 


129 Sn(f + ) 


5.57 


5.9(3) 


49.8 


34(12) 




129 Sn(i + ) 


5.80 


5.5(1) 


38.1 


54(6) 


2 ° 5 Hg(r) 


205 T1(I + ) 


5.37 


5.257(11) 


62.3 


71.3(9) 




205 rpl/ 1 +\ 

Ll \22 ' 


6.77 


7.03(25) 


12.5 


9(3) 




205 T1(§ + ) 


7.26 


6.51(21) 


18.9 


17(4) 




205 T1 / 3 +\ 
Ll \22 1 


6.32 


7.61(22) 


10.3 


5(1) 




205 T1(| + ) 


8.16 


8.70(21) 


1.91 


1.3(3) 


206 Hg(0+) 


206 T1(0") 


5.42 


5.41(6) 


59.2 


60(4) 




206 T1(I-) 


5.18 


5.24(10) 


77.6 


73(8) 




206 T1(1 2 -) 


5.68 


5.67 


43.6 


44.3 


207 T1(± + ) 


207 Pb(I") 


5.14 


5.108(6) 


81.7 


84.5(6) 




207 Pb(D 


6.18 


6.157(22) 


24.7 


25.3(6) 


206 T1(0-) 


206 Pb(0+) 


5.42 


5.1775(13) 


52.4 


78.0(1) 




206 Pb(0+) 


5.18 


5.99(6) 


32.4 


31(2) 




206 Pb(2+) 


5.68 


8.60(3) 


1.87 


1.52(5) 


205 Au(f + ) 


205 Hg(i-) 


6.79 


5.79(9) 


12.1 


39(4) 




205 Hg(f:) 


7.33 


6.43(11) 


6.5 


18(2) 




205 Hg(f-) 


5.82 


6.37(12) 


37.3 


20(3) 



For the nuclei in the vicinity of N = 126 we find, how- 
ever, a noticeably larger scatter between calculation and 
data. Satisfyingly there are no systematic deviations. 
With the exception of the two 205 Au decays, where our 
calculation overestimates (to the 5/2~ state in 205 Hg) 
or underestimates (to the l/2~ state) the average shape 
factor roughly a factor 9, we generally find agreement of 
our calculated C(W) with data within a factor of 4. As 
already observed while determining the quenching factor 
for GT transitions in shell model calculations [54| the de- 
scription of a decay between specific states is noticeably 
more sensitive to nuclear structure effects than global 
quantities like half-lives or total strengths. Hence we ex- 
pect that our prescription of quenching for first-forbidden 
transitions yields a fair description of the N = 126 half- 
lives. 

As already stressed by Warburton [5(| the relativistic 
matrix element (Eq. 8h) is enhanced compared to the 
other first-forbidden transitions. We confirm this finding 
as the value for g(£'i>) is noticeably larger than the other 
quenching factors. 

Having determined the quenching of first-forbidden 
transitions, we adjust the quenching of the Gamow- Teller 
transition to the half-life of 130 Cd which is expected 
to decay dominantly by Gamow- Teller. This is indeed 
borne out in our calculation. Using the quenching fac- 



tor qcT = 0.66 we reproduce the measured half-life using 
both GT and first-forbidden transitions. The latter con- 
tribute about 13% to the half-life and hence are small, 
but not negligible. Our factor qoT is only slightly smaller 
than the customary quenching value of 0.7. 

All half-lives presented in the next sections for the r- 
process waiting point nuclei have been obtained using the 
quenching factors for Gamow- Teller and first-forbidden 
transitions derived above. 

Before we present our results we have to discuss an- 
other potential shortcoming of our calculations of first- 
forbidden transitions and how we will handle it. A com- 
pletely converged calculation of the first-forbidden tran- 
sition strength in our chosen model space is prohibited 
due to computational limitations. We have derived the 
strength within the Lanczos scheme using 100 iterations. 
As a consequence the lowest states are converged and 
correspond to physical states, while the Lanczos states 
at higher excitation energies are unphysical and repre- 
sent strength per energy interval. Furthermore, their en- 
ergy positions depend on the sum rule (pivot) state used 
for the calculation of the strength function. For exam- 
ple, starting from the pivot state of the x operator one 
obtains different energy positions for the non-converged 
states than starting from the sum rule state of the u 
operator. As we need to compute superpositions of oper- 
ators like u + x, where both the magnitude and the phase 
of the individual operators matter, we have followed the 
same procedure as is used in shell-model calculations of 
double-beta decays [32j . Hence we start with an arbitrary 
sum rule state that can be any linear combination of op- 
erators of the same rank. During the Lanczos iteration 
procedure we compute the overlaps with the individual 
operators. This iteration procedure is stopped when at 
least 80% of the total strength for each individual oper- 
ator is recovered. 

To illustrate this point, we have performed calcula- 
tions for 199 Ta in the model space defined above and have 
used two different linear combinations of rank 1 opera- 
tors, £'y — £(x + u) and (x + u), as pivots for the Lanczos 
calculations of the operators x and u. For the combina- 
tion £'y — l;(x + u) we obtain 92% and 80% of the total 
strength for the operators x and u, respectively, while for 
the combination x + u we recover 85% and 95% after 100 
iterations. We stress that the effect which these short- 
comings have on the first-forbidden half-lives is mildened 
as the contributions of the low-lying states, which are 
converged in our Lanczos scheme, are strongly enhanced 
by the phase space energy dependence. 

To quantify the potential uncertainty in our first- 
forbidden half-lives, we have performed again calcula- 
tions for 199 Ta in the model space as defined above, how- 
ever, allowing only one proton pair to be broken in our 
seniority scheme. This truncated space allows for the 
calculation of a fully converged first-forbidden strength 
distribution in the Qp window. Fig. [5] compares the 
partial decay rates to the various states in the daugh- 
ter nucleus obtained in the fully converged calculation 
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FIG. 2. Partial decay rates for Ta calculated from first- 
forbidden rank 1 transitions only. The top (a) and medium 
(b) panel shows the decay rates obtained by using the lin- 
ear combinations £'y — £(a; + it) and (x + u), respectively, as 
pivot elements to calculate the contributions of rank 1 opera- 
tors within a Lanczos scheme with 100 iterations. The bottom 
panel (c) shows the partial decay rate using a Lanczos scheme 
with 300 iterations in which the states within the Qp window 
are converged. These shell model calculations have been per- 
formed in a truncated model space compared to other studies 
of the N — 126 isotones. 



(with 300 Lanczos iterations) with those where the rank 
1 contributions were derived using the linear combina- 
tions £'y — £(cc + u) and (x + u) as pivots for the Lanc- 
zos scheme with 100 iterations. As expected, the lowest 
Lanczos states are converged and hence the calculated 
strength is the same for either choice of pivot combina- 
tions. The Lanczos scheme with 100 iterations is not suf- 
ficient to converge the states for excitation energies larger 
than about 2.5 MeV. They represent unphysical states as 
discussed above. Obviously further iterations lead to a 
stronger fragmentation of the strength in a small energy 
interval around the unphysical states. The energy inter- 
val is small enough that this redistribution of strength 
due to different phase space weighting has negligible ef- 
fect on the half-life. But also the interference contribu- 
tions lead to rather mild differences between the trun- 
cated and the converged calculations. We find a partial 
half-life due to rank 1 first-forbidden operators of 655 ms 



and 716 ms when calculating the rank 1 operators from 
the linear combinations £'y — £,(x + u) and (x+u), respec- 
tively, while the partial half-life in the converged study is 
651 ms agreeing within 10 percent with the two approx- 
imate calculations. In the following we will calculate the 
contributions from the first-forbidden rank 1 operators 
using a Lanczos scheme with the pivot state — £(x + u) 
and 100 iterations. 



IV. HALF-LIVES OF THE N = 50 WAITING 
POINT NUCLEI 



> 




a io 




26 28 30 

Charge number Z 

FIG. 3. (Color online) Comparison of Qp values of the N = 
50 isotones between experimental data [55| and theoretical 
results from FRDM [lj] and the present shell model approach. 

To calculate half-lives, a good description of the tran- 
sition matrix elements and also of the Qp values is re- 
quired. As is demonstrated in Fig. [3] our shell model cal- 
culation reproduces the Qp-vahies as given in the Audi- 
Wapstra compilation well (55[. Hence we will use the 
shell model Qp values in the following calculation of the 
half-lives and /3-delayed neutron emission probabilities 
for the N = 50 waiting point nuclei. 

As is shown in Table ILT1 and in Fig. [4] the shell model 
half-lives agree quite well with the data, although they 
overestimate the ones of 82 Ge and 79 Cu by about 50%. 
Nevertheless the agreement is significantly better than 
obtained based on the global FRDM and ETFSI models. 
The HFB results [l6j], which are restricted to the decay 
of even-even nuclei, are very similar to the shell model 
results, except for the half-life of the double-magic nu- 
cleus 78 Ni. Here only the shell model reproduces the 
measured value [7|, while all other models predict a sig- 
nificantly longer half-life. This underlines again the fact 
that many-body configuration mixing is needed to repro- 
duce the cross-gap correlations in double-magic nuclei. 
Similar results have been found in studies of the isotope 
shifts in calcium |57l| or the Ml strength distributions in 
argon isotopes [58] . 
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TABLE II. Comparison of the present shell model half-lives 
and the ones of reference I21JI. All half-lives are in s. 



Nucleus 




Half-Life (s) 




Expt. 


present shell model (ref [2111 



2 Ge 



82 ( 

81 Ga 
80 Zn 
79 Cu 
78 Ni 
77 Co 
76 Fe 
75 Mn 
74 Cr 



4.55 ±0.05 

1.217 ±0.005 

0.545 ±0.016 

0.188 ±0.025 

11+ 01 
u - 11 -o.o6 



6.90 
1.03 
0.53 
0.27 
0.15 
0.016 
8.26 x 10" 
3.66 x 10" 
2.23 x 10" 



2.057 
0.577 
0.432 
0.222 
0.127 
0.016 
7.82 x 10" 
3.52 x 10" 
2.07 x 10" 




26 28 30 

Charge number Z 

FIG. 4. (Color online) Comparison of half-lives of the N = 50 
isotones between experimental data and theoretical results 
from FRDM+QRPA 14], HFB±QRPA QH, DF3+QRPA 
20, 56] and the present shell model approach. 



The contribution of first-forbidden transitions to the 
N = 50 half-lives is shown in Fig. [5J For the decay of 
the nuclei with Z > 28 the probability is very small (less 
than 5%). However, first-forbidden transitions contribute 
about 25% to the 77 Co decay, while they are smaller, 
but still sizable for the decay of the nuclei with charge 
numbers Z — 24-26. 

To understand this behavior, we note that first- 
forbidden contributions are related to the transition from 
a (79/2 neutron orbital to a f?/2 proton orbital for the 
rank 1 operators and to a f$/2 proton orbital for rank 
2 operators. (There are no contributions from rank 
operators in our model space.) In the simple Indepen- 
dent Particle Model the level gets completely occu- 
pied for 78 Ni and consequently this transition is Pauli- 
blocked for N = 50 nuclei with Z > 28. In the shell 
model, the blocking is partially removed by configura- 
tion mixing, but the importance of the first-forbidden 




26 28 30 

Charge Number Z 

FIG. 5. (Color online) Percentage of the contributions from 
first-forbidden transitions to the half-lives of the N= 50 iso- 
tones from FRDM±QRPA 14], DF3+QRPA [UIH and the 
present shell model. 



transitions stay low. For these nuclei they are nearly ex- 
clusively due to contributions from the rank 2 operators. 
For the nuclei with Z < 28 the proton / 7 / 2 orbital is 
not fully occupied and first-forbidden transitions due to 
the rank 1 operators are possible. They are noticeably 
larger those of the rank 2 operators. Hence the total first- 
forbidden strength is significantly larger for nuclei with 
Z < 28 than for the nuclei with Z > 28. Furthermore it 
increases with decreasing charge number due to the de- 
population of the proton f^M orbital in the daughter nu- 
cleus. However, also Gamow- Teller transitions from the 
neutron / 5 / 2 orbital into the proton / 7 / 2 orbital become 
unblocked. Hence reducing the charge number, increases 
both the GT and first-forbidden transitions due to de- 
creasing Pauli-blocking of the dominant transitions into 
the ffj2 orbital. However, the relative decrease for the 
GT half-lives with decreasing charge number is stronger 
than for the first-forbidden transitions. This is related to 
phase space. Examples of differential decay rates as func- 
tion of excitation energies for nuclei 74 Cr, 77 Co, 78 Ni and 
82 Ge are shown in Figs. [6] From these figures we note 
that the first-forbidden transitions are dominantly pro- 
ceeding to states in the daughter at low excitation ener- 
gies (usually up to 2-2.5 MeV) for the N = 50 nuclei with 
Z < 28, while the GT transitions go to states with ex- 
citation energies of order 5-7 MeV, simply reflecting the 
fact that it is energetically more favorable to have a f 5 / 2 
neutron hole and a closed 179/2 shell, than having a hole 
in the g$/2 orbital. As the energy gain in the transitions 
is smaller for the GT transitions, they are more sensitive 
to the increase of the Qp value with decreasing charge 
number. This explains why the relative contribution of 
first-forbidden transitions decreases with reduced charge 

Ni the 



78 



number below the double-magic Ni. Ab 
GT transitions proceed to daughter states at relatively 
low excitation energies. (Fig. [B] shows the differential 
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FIG. 6. (Color online) Partial decay rates including GT and FF transitions for the N = 50 isotones 74 Cr (a), 78 Ni (b), 77 Co 
(c), and 82 Ge (d). 




decay rates for 82 Ge as an example.) As a consequence 
first forbidden transitions, due to their smaller transition 
matrix elements, cannot compete with GT transitions. 



We have used the shell model neutron separation en- 
ergies to calculate the /3-delayed neutron emission prob- 
abilities, which are shown in Fig. [7] As demonstrated 
above GT transitions dominate the decays of the N = 50 
r-process waiting points. As for the nuclei with Z < 29 
these transitions connect mainly to daughter states at 
excitation energies above the respective neutron thresh- 
old, the decay is accompanied by neutron emission with 
a high probability. The probability comes close to 100% 
for the decays of 55 Mn and 54 Cr associated with the very 
low neutron separation energies in the daughters. For the 
nuclei with Z > 29, the GT transitions reside at lower 
excitation energies in the respective daughters, while the 
neutron thresholds increases. As a result the neutron 
emission probability is strongly reduced in these nuclei. 



V. HALF-LIVES OF THE N = 82 WAITING 
POINT NUCLEI 

The present interaction and model space, based on a 
88 Sr core, are not the same as used in Ref. [l9[ to cal- 
culate Qp values and Gamow- Teller strength functions. 
However, we stress that the present shell model calcula- 
tion gives very similar results to the ones of Ref. [Tj| . In 
particular we reproduce the experimentally available Qp 
values very well, as is shown in Fig. [3] This figure also 
shows that the agreement of the Qp values obtained in 
other models is usually not as good as by the shell model 
results. We will in the following use the shell model Qp 
values for the calculation of the half-lives. 

The 1/2 - isomer in 131 In corresponds approximately 
to a 132 Sn configuration with a hole in the lpi/2 orbital. 
Our calculation reproduces the energy of the isomer at 
0.302 MeV. (This quantity was one of the experimental 
ingredients to which the interaction has been adjusted.) 

The calculated half-lives for the N = 82 waiting point 
nuclei are summarized in Table Mil and are compared 
to data and to previous theoretical estimates in Fig. 
[5] Compared to experiment, the half-life of 131 In is 
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FIG. 7. (Color online) Neutron emission probability for N= 
50 isotones from FRDM+QRPA [3, DF3+QRPA [U [U 
the present shell model and experiment [53| . 
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FIG. 9. (Color online) Comparison of half-lives of the N = 82 
isotones as calculated in the FRDM, HFB, DF3+QRPA [2(| 
and the present shell model approaches with data. 
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FIG. 8. (Color online) Comparison oiQp values of the N = 82 
isotones between theoretical results from the FRDM [hH ]. 
HFB [l6l l , th e present shell model approaches and experimen- 
tal data \M- 



well reproduced, while the one for 129 Ag is somewhat 
too long. This shortcoming had already been observed 
in the previous shell model calculations. In fact, the 
present shell model results, including contributions from 
first-forbidden transitions, agree very well with the shell 
model results of Ref. 19]. This, however, does not 
mean that first-forbidden transitions are negligible. As 
is shown in Fig. [10] first- forbidden transitions contribute 
about 13% to the half-life. However, this value is nearly 
the same for all N = 82 waiting point nuclei explaining 
the similarity between the present shell model results to 
those of Ref. Only for 131 In, first-forbidden tran- 

sitions contribute somewhat more, resulting in a slightly 
smaller half-life than in the shell model study based 
solely on Gamow- Teller transitions. We note that our 



TABLE III. Comparison of the present shell model half-lives 
and the ones of reference pj| with experiment [Hi], [5j| . All 
half-lives are in ms. 



Nucleus 




Half-Life (ms) 


Expt. 


present 


shell model (ref [19]) 


131 In 


280 ± 30 


247.53 


260 


130 Cd 


162 ± 7 


164.29 


162 


129 Ag 




69.81 


70 


128p d 




47.25 


46 


127 Rh 




27.98 


27.65 


126 Ru 




20.33 


19.76 


125 T( , 




9.52 


9.44 


124 Mo 




6.21 


6.13 



prediction of an 18% contribution stemming from first- 
forbidden transitions to the decay of the 131 In ground 
state is in agreement with the experimental limit of 
< 20%. We further add that we calculate a Gamow- 
Teller contribution to the half-life of the 1/2" isomer 
in 131 In which is less than 1%, confirming our assump- 
tion to fix the quenching of the first-forbidden transition 
to this decay. Fig. [IT] shows the partial decay rates 
to different final states for the nuclei 124 Mo, 126 Ru and 
128 Pd. We note that Gamow- Teller transitions are larger 
than first-forbidden transitions, which, however, proceed 
to levels at lower excitation energies which enhances them 
by phase space. 

The /3-delayed neutron emission probabilities, i.e. the 
probabilities that the decay leads to states in the daugh- 
ter nucleus above the neutron separation threshold and 
hence is followed by the emission of a neutron, is obvi- 
ously sensitive to a good description of both the neutron 
separation energies and the /3 strength functions in the 
Qp window. As has been stressed in Refs. [13, Qjl the 
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FIG. 10. (Color online) Percentage of the contributions from 
first-forbidden transitions to the half-lives of the N= 82 iso- 
tones from FRDM+QRPA 14], DF3+QRPA [H[5(| and the 
present shell model data. 



improved description of correlations in shell model calcu- 
lations gives a more realistic account of the fragmenta- 
tion of the strength function than is obtained in QRPA 
studies. Fig. [T3] compares the present shell model prob- 
abilities to those obtained in Ref. We find that the 
inclusion of first-forbidden transitions leads only to minor 
changes. A detailed comparison of the shell model results 
[l^ | to those obtained in other theoretical approaches is 
given in [l7l[l9j. 



VI. HALF-LIVES OF THE N = 126 WAITING 
POINT NUCLEI 

Fig. [13] compares the calculated Qp values of the 
N=126 isotones with other theoretical models. While 
the general trend of the Qp is quite similar than obtained 
in the FRDM model, the shell model values are slightly 
smaller than those from the FRDM model. As, however, 
no experimental data exist for these very neutron-rich 
N = 126 nuclei, it is not possible to decide which Qp are 
more realistic. In the following we will use the shell model 
values to calculate the half-lives and /3-delayed neutron 
emission probabilities for the N — 126 r-process waiting 
point nuclei. 

The shell model half-lives are listed in Table ITVl and are 
compared to other theoretical predictions in Fig. Q3J Al- 
though recently researchers at GSI have been successful 
to measure half-lives of nuclei close to N = 126 with 
charge numbers below lead supplying important 

constraints about the half-life trend towards the r-process 
nuclei, experimental data for the N — 126 r-process nu- 
clei do yet not exist. Hence our results can only be com- 
pared to other theoretical predictions. We note that the 
present half-lives for Z > 70 are faster, by about a factor 
of two, than those obtained by Borzov within an QRPA 



TABLE IV. Comparison of the present shell model half-lives 
and the one of reference [601. All half-lives are in ms. 



Nucleus 



Half-Life (ms) 



present 



SM (ref [60J 



:J Ta 



199 n 

198 Hf 

197 Lu 

196 Yb 

195 Tm 

194 Er 
193 



Ho 

2 Dy 



286.17 
193.28 
107.85 
68.98 
36.03 
24.58 
13.58 
10.10 



278.88 
129.65 
84.81 
44.18 
29.49 
18.11 
10.94 
7.75 



approach on top of the density functional DF3, sho wing , 
however, a similar dependence with charge number |56| . 
Adopting a different parameterization, Borzov has also 
calculated half- lives for N = 126 isotones with Z < 70 
which are slightly faster than the shell model values [62j ■ 
The shell model half-lives are noticeably faster than those 
predicted previously by global models, e.g. by the QRPA 
calculation on top of the microscopic-macroscopic FRDM 
or ETFSI approaches. 

Recently Suzuki et al. [60] have presented the first shell 
model half-lives for N — 126 r-process nuclei, includ- 
ing both GT and first-forbidden contributions. However, 
our present model space including the (0<?7/ 2 , 1^5/2,3/21 
O/in/2, 2s 1 / 2 ) proton orbits is noticeably larger than the 
one used in Ref. [6(| (the ld 3 / 2 , 0/i n / 2 , 2s!/ 2 proton or- 
bits). Relatedly the two shell model calculations differ in 
the residual interaction and additionally in the adopted 
quenching scheme for first-forbidden transitions. Never- 
theless, as is shown in Table IIVI and Fig. [TH both shell 
model calculations predict very similar half-lives for the 
N = 126 nuclei. Both studies do not predict the strong 
odd-even staggering in the half-lives as observed in the 
QRPA results on top of the FRDM model. 

As already noted in Ref. |20| based on the density 
functional calculations, first-forbidden transitions are ex- 
pected to contribute significantly to the half-lives of the 
N = 126 r-process nuclei. This finding is supported by 
our shell model calculations (Fig. fT5|) . One observes an 
increasing contribution from the first-forbidden transi- 
tions with increasing proton number. In fact, for nu- 
clei with proton number Z > 70, contributions from 
first-forbidden transitions to the half-life dominate over 
Gamow- Teller transitions. This behavior can be under- 
stood by inspecting the partial decay rates arising from 
Gamow- Teller and first-forbidden transitions which are 
shown in Fig. [16] for selected nuclei. We note that 
Gamow- Teller transitions are related to the change of 
a neutron in the Oft.9/2 orbit to a proton in the 0/in/ 2 
orbit which, however, is fragmented over several states 
in the daughter nucleus due to correlations. Neverthe- 
less, for the nuclei studied here these final proton states 
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FIG. 11. (Color online) Partial decay rates including GT and FF transitions for the N = 82 isotones 124 Mo (a), 126 Ru (b) and 
128 Pd (c). 
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FIG. 12. (Color online) /3-delayed neutron emission probabil- 
ity for selected N = 82 r-process nuclei from FRDM+QRPA 
pj], DF3+QRPA [HE!], the present shell model and exper- 
iment 15311. 
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FIG. 13. (Color online) Comparison of Qr values of the N = 
126 isotones as calculated in the FRDM [lj] and the present 
shell model approaches. 



reside at moderately high excitation energies around 3 
MeV, while first-forbidden transitions connect to excited 
states at lower excitation energies. With increasing pro- 
ton number, more protons occupy the final 0/in/ 2 orbit 
and the GT transitions get gradually Pauli blocked. This 
explains why the GT strength gets strongly reduced with 
increasing proton number. Actually also first-forbidden 
transitions get blocked with increasing proton number, 
which is, however, a significantly milder effect as for the 
GT transitions. We note that for all nuclei studied here 
first-forbidden transitions are mainly mediated by rank 
and 1 operators (the latter contributes about 70% to the 



forbidden strength), while the contribution arising from 
rank 2 operators are very small. Due to its larger sen- 
sitivity to phase space (~ Q ) the relative contribution 
of rank 2 transitions increases slightly with decreasing 
charge number. For the N — 126 nuclei studied here, the 
depopulation of the 0hn/ 2 proton orbital with decreasing 
charge number, which increases both the Gamow- Teller 
and first-forbidden transitions, dominates the trend ob- 
served in our half-lives. Changes in phase space and pair- 
ing have lesser effects on our half-life systematics. 

Gamow- Teller transitions proceed to final states 
mainly above the neutron threshold and hence are ac- 
companied by neutron emission, while the final states 
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FIG. 14. (Color online) Comparison of half-lives of the 
N = 126 isotones as calculated in the FRDM+QRPA, 
DF3+QRPA(I) [U, DF3+QRPA(II) [H| and the present 
shell model approaches [ool ] . 
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FIG. 15. (Color online) Percentage of the contributions from 
first-forbidden transitions to the half-lives of the N= 126 iso- 
tones are compared with results from DF3+QRPA(I) [20l |. 
DF3+QRPA(II) 61] and shell model approaches [U- 



populated by first-forbidden transition predominantly re- 
side below the neutron threshold. Hence we expect from 



the Z-dependence of the GT and first-forbidden tran- 
sitions that the /3-delayed neutron emission probability 
decreases with increasing proton number. This is indeed 
confirmed by Fig. 1171 The striking odd-even staggering 
is related to pairing which reduces the neutron threshold 
energies in odd-odd daughter nuclei relatively to odd-A 
nuclei, but does basically not affect the strength distri- 
butions, as discussed in [63|. The QRPA/FRDM neu- 
tron emission probabilities show a rather abrupt increase 
at Z — 68 which is likely due to the fact that QRPA 
calculations show significantly less fragmentation of the 
strength than shell model studies and that, for Z < 69, 
the few dominant transitions reside above the neutron 
threshold. 



VII. SUMMARY AND CONCLUSIONS 

We have calculated the half-lives and /3-delayed neu- 
tron emission probabilities of the r-process waiting point 
nuclei with magic neutron numbers N — 50, 82, and 
126 within the framework of the large-scale shell-model. 
The calculations include contributions both from allowed 
Gamow- Teller and first-forbidden transitions. We find 
good agreement with the existing experimental data: i.e. 
the half-lives for the TV = 50 nuclei with charge num- 
bers Z = 28-32 and for the N = 82 nuclei 129 Ag and 
130 Cd. In our calculations first-forbidden transitions sig- 
nificantly reduce the half-lives of the N = 126 waiting 
point nuclei, while they have a smaller effect on the half- 
lives of the N = 50 and 82 r-process nuclei. 
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